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EXECUTIVE  SUMMARY 


A.  OBJECTIVE 

The  objective  of  this  study  was  to  gain  some  insight  into  failure  mechan¬ 
isms  of  concrete  at  strain  rates  associated  with  high  intensity  loadings  from 
conventional  explosives.  To  this  end,  a  comprehensive  numerical  analysis  of 

splitting-tensile  and  direct  tension  tests  of  plain  concrete,  performed  at 

2 

strain  rates  between  1  and  10  per  second  on  a  Split  Hopkinson  Pressure  Bar 
(SHPB),  was  conducted. 

B.  BACKGROUND 

The  understanding  of  material  response  to  high  amplitude,  short-duration, 
impulse  loads  generated  in  a  weapons  environment  is  an  important  problem  in 
protective  construction  design  and  analysis.  To  model  the  response  in  the 
laboratory  requires  that  the  environment  must  reflect  the  type  of  confinement, 
magnitude  of  stress  change,  and  the  time  scale  of  loading  anticipated  in  the 
field.  The  Split  Hopkinson  Pressure  Bar  (SHPB)  technique  can  produce  the 
required  environments  in  the  laboratory. 

However,  several  significant  shortcomings  are  associated  with  SHPB 
experiments.  First,  it  is  not  possible  to  accurately  determine  the  stress 
condition  in  the  specimen  at  failure  from  the  available  data,  and  second,  it 
is  frequently  not  possible  to  ascertain  the  mode  of  failure  in  the  specimen. 
Therefore,  a  comprehensive  numerical  analysis  was  conducted  on  various  SHPB 
experiments  to  gain  some  insight  into  those  phenomena. 

C.  SCOPE 

Two  different  types  of  SHPB  experiments  were  simulated  in  the  numerical 
analyses:  (1)  splitting-tensile  tests  and  (2)  direct  tension  tests.  In  the 
splitting-tensile  analyses,  the  numerical  model  included  the  2-inch  (51  rrni) 
diameter  cylindrical  specimen  and  a  10-inch  (25.4  cm)  length  of  the  trans¬ 
mitter  bar.  Three  different  load  cases  were  investigated.  In  the  direct 
tension  study,  both  a  square  notch  and  saddle  notch  specimen  were  analyzed. 

The  numerical  model  included  the  entire  lengths  of  the  incident  and  trans¬ 
mitter  bars,  in  addition  to  the  2-inch  (51  mm)  diameter  specimen.  One  load 
case  for  each  specimen  type  was  investigated. 
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D.  METHODOLOGY 

Because  of  the  dynamic  nature  of  the  impulse  loading  associated  with  the 
SHPB  experiments,  and  the  highly  nonlinear  behavior  of  the  concrete  test 
specimens,  the  finite  element  method  (FEM)  of  analysis  was  employed  in  the 
research  effort  through  implementation  of  the  ADINA  computer  programs.  The 
analyses  were  conducted  on  two-dimensional,  axisymmetric  models  comprised  of 
nine-node  isoparametric  finite  elements.  Both  linear  and  nonlinear  analyses 
were  performed. 

E.  TEST  DESCRIPTION 

The  high  strain  rate  SHPB  experiments  were  performed  using  the  Air  Force 
Engineering  and  Services  Center  (AFESC)  50.8  mm  diameter  SHPB  by  AFESC  per¬ 
sonnel.  All  numerical  analyses  were  conducted  on  the  Alabama  Supercomputer 
Network  (ASN)  Cray  X-MP/24  Supercomputer  by  Auburn  University  personnel. 

F.  RESULTS 

Both  an  elastic  and  an  inelastic  concrete  material  model  were  employed 
in  all  numerical  analyses.  For  the  splitting-tensile  study,  the  results  of 
the  linear  analyses  indicate  that  the  dynamic  stress  distribution  in  the 
cylinder  behind  the  initial  stress  wave  is  identical  to  that  exhibited  in 
static  analyses.  The  linear  results  also  indicate  that  the  maximum  tensile 
stress  always  occurs  at  the  center  of  the  cylinder.  This  observation  is  like¬ 
wise  consistent  with  the  results  of  the  static  analysis. 

The  mode  of  failure  predicted  by  the  nonlinear  analysis  differs  from 
that  suggested  by  the  results  of  the  linear  analysis.  In  all  three  load  cases, 
the  initiation  of  first  cracking  is  not  at  the  center  of  the  cylinder  (as  the 
results  of  the  linear  analysis  indicate),  but  at  an  approximate  distance  of 
0.2D  from  the  top  of  the  cylinder.  The  sequence  of  failure  for  all  load  cases 
is  essentially  the  same:  initiation  of  first  cracking  at  location  0.2D  from 
the  top,  and  subsequent  propagation  of  the  cracks  in  both  directions  along  the 
vertical  centerline  of  the  cylinder  toward  the  top  and  bottom  surfaces.  Some 
mid-diameter  crack  bifurcation  occurs  as  the  load  rate  is  increased.  This 
prediction  has  been  verified  by  observation  of  experimental  results  using 
high  speed  photography. 
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For  the  direct  tension  study,  the  results  of  the  linear  analyses  indicate 
the  development  of  high  stress  concentrations  at  the  root  of  the  notch  in  the 
square-notch  specimens,  and  at  the  apex  of  the  notch  in  the  saddle-notch 
specimens.  These  results  suggest  that  first  cracking  will  begin  at  these 
locations  of  high  stress  concentration  factors  and  that  failure  will  occur  on 
vertical  planes  passing  through  these  locations. 

The  mode  of  failure  predicted  by  the  nonlinear  analysis  differs  from 
that  suggested  by  the  results  of  the  linear  analysis.  In  the  square-notch 
specimen,  first  cracking  occurs  at  the  root  of  the  notch.  This  is  consistent 
with  the  stress  concentration  predictions  from  the  linear  analysis.  However, 
eventual  failure  of  the  specimen  occurs  on  a  vertical  plane  adjacent  to  the 
face  of  the  incident  bar.  In  the  saddle-notch  specimen,  first  cracking  occurs 
at  a  transverse  section  in  the  specin«n  next  to  the  indicent  bar.  Almost 
simultaneously,  cracks  develop  in  the  apex  of  the  notch.  Eventual  failure 
is  along  the  transverse  section  adjacent  to  the  incident  bar. 

G.  CONCLUSIONS 

In  the  case  of  the  splitting-tensile  tests,  it  can  be  concluded  that  the 
nature  of  the  failure  mode  is  directly  affected  by  the  rate  of  loading.  For 
a  relatively  low  load  rate,  the  failure  mode  manifests  itself  as  a  single 
crack  propagating  along  the  vertical  centerline  of  the  cylinder.  However,  for 
increasingly  higher  load  rates,  the  mode  of  failure  is  characterized  by 
several  bifurcations  in  the  primary  crack  pattern.  The  higher  the  load  rate, 
the  more  pronounced  are  the  bifurcations. 

In  the  case  of  the  direct- tension  tests,  the  results  of  the  linear  analyses 
indicate  high  stress  concentration  factors  in  the  vicinity  of  the  notches. 

These  results  suggest  a  failure  in  a  transverse  plane  passing  through  the 
notches.  The  nonlinear  analyses,  however,  predict  failure  in  a  transverse 
plane  near  the  end  of  the  specimen  next  to  the  incident  bar.  The  reason  for 
this  is  that  the  load  rate  is  so  high,  that  the  tensile  limit  of  the  material 
is  reached  at  the  end  of  the  specimen  (adjacent  to  the  incident  bar)  before 
any  significant  stresses  can  develop  on  the  transverse  plane  passing  through 
the  notch. 
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H.  RECOMMENDATIONS 

The  failures  predicted  in  both  the  splitting-tensile  tests  and  the  direct 
tension  tests  are  highly  sensitive  to  the  rate  of  loading.  Therefore  it  is 
recommended  that  additional  analyses  be  conducted  at  a  wide  range  of  load 
rates  to  quantify  the  relationship  of  load  rate  to  mode  of  failure.  It  is 
also  apparent  from  the  results  of  the  analyses  that  material  strain  rate 
effects  will  delay  the  time  of  failure,  allowing  the  specimen  to  be  subjected 
to  a  higher  load,  thus  possibly  affecting  the  failure  mode.  Therefore,  it  is 
further  recommended  that  additional  numerical  analyses  be  conducted  to  investi¬ 
gate  material  strain  rate  effects  on  the  mode  of  failure. 

Finally,  the  notches  in  the  direct  tension  specimens  analyzed  in  this 
study  were  relatively  shallow.  It  is  recommended  that  specimens  with  deeper 
notches  be  analyzed,  both  experimentally  and  numerically,  to  quantify  the 
effect  of  notch  depth  on  the  mode  of  failure. 
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SECTION  I 


INTRODUCTION 

A.  OBJECTIVE 

A  comprehensive  numerical  analysis  of  splitting-tensile  and  direct 

2 

tension  tests  of  plain  concrete,  performed  at  strain  rates  between  1  and  10 
on  a  Split  Hopkinson  Pressure  Bar  (SHPB),  was  conducted  to  ascertain  the 
states  of  stress  in  the  concrete  specimens  at  failure  and  to  identify  the 
modes  of  failure. 

B.  BACKGROUND 

1.  Split  Hopkinson  Pressure  Bar  (SHPB) 

The  understanding  of  material  response  to  high-amplitude,  short- 
duration,  impulse  loads  generated  in  a  weapons  environment  is  an  important 
problem  in  protective  construction  design  and  analysis  (Reference  1).  To 
model  the  response  in  the  laboratory  requires  that  the  environment  must  re¬ 
flect  the  type  of  confinement,  magnitude  of  stress  change,  and  the  time  scale 
of  loading  anticipated  in  the  field  (Reference  2).  The  Split  Hopkinson 
Pressure  Bar  (SHPB)  technique  (Reference  3)  can  produce  the  required  environ¬ 
ments  in  the  laboratory. 

During  the  past  several  years  researchers  have  demonstrated  that  the 
SHPB  technique  can  determine  the  dynamic,  high  stress  and  strain  rate  of  soil 
(References  4,  5,  6,  7,  and  8)  and  concrete  (References  9,  10,  11,  and  12). 
Although  the  conditions  of  the  experiment  have  been  restrictive  (e.g.,  condi¬ 
tions  of  uniaxial  strain),  this  technique  has  significantly  extended  the 
stress  and  strain-rate  regimes  over  which  dynamic  material  properties  can  be 
investigated. 

Hopkinson  (Reference  13)  introduced  the  concept  of  using  a  cylindri¬ 
cal  bar  for  evaluating  material  response  to  impulse  loads.  The  apparatus  con¬ 
sisted  of  a  long  cylindrical  bar  with  a  time  piece  of  the  same  diameter  and 
material  attached  by  magnetic  attraction  to  one  end.  By  propagating  a  com¬ 
pressive  wave  down  the  bar  and  capturing  the  momentum  transferred  to  time 
pieces  of  different  lengths,  an  approximate  stress-ti.ne  curve  could  be  con¬ 
structed.  Davies  (Reference  14)  improved  the  experimental  technique  by  intro¬ 
ducing  electrical  condenser  units  to  measure  the  displacement  at  the  surface  of 
the  bar  caused  by  the  propagating  wave.  In  addition,  Davies  developed  a 
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theoretical  foundation  of  dispersion  phenomenon  and  established  the  accuracy  of 
the  technique  when  assuming  one-dimensional  wave  propagation  in  the  bar.  Using 
this  framework,  Kolsky  (Reference  15)  modified  the  technique  to  permit  the 
dynamic  response  of  a  material  to  be  measured  indirectly  by  placing  a  specimen 
between  two  bars  fitted  with  condenser  microphones  for  data  recording.  Assum¬ 
ing  that  a  uniform  distribution  of  stress  existed  along  the  longitudinal  axis 
of  the  specimen,  Kolsky  developed  relationships  to  compute  the  average  stress, 
strain,  and  strain-rate  response  of  the  specimen.  This  technique  is  now  known 
as  the  Kolsky  technique  or  Split  Hopkinson  Pressure  Bar  technique. 

An  illustration  of  the  SHPB  device  is  shown  in  Figure  1.  The  device 
is  operated  by  the  Engineering  and  Services  Laboratory,  Air  Force  Engineering 
and  Services  Center,  Tyndall  AFB,  Florida.  The  pressure  bars  are  constructed 
of  PH  13-8  MO  stainless  steel.  Each  pressure  bar  is  2.0  inches  (51  mm)  in 
diameter.  The  lengths  of  the  incident  and  striker  bars  are  12  and  11  feet 
(3.66  and  3.35  m),  respectively.  Striker  bar  lengths  of  4,  6,  and  8  inches 
(102,  153,  and  203  mm),  are  available.  The  loading  compressive  stress  wave  is 
initiated  by  the  impact  of  the  striker  bar  (which  is  propelled  by  the  gas  gun) 
on  the  incident  bar  (Figure  2).  The  amplitude  of  the  incident  stress  pulse  is 
determined  by  the  impact  velocity  and  material  properties  of  the  striker  bar, 
while  the  duration  of  the  pulse  depends  on  the  length  and  wave  speed  of  the 
striker  bar  (Reference  16). 

The  incident  stress  wave  (oj)  generated  in  the  incident  bar  travels 
down  the  bar  and  is  recorded  at  Strain  Gage  A  (Figure  2),  is  partially  reflected 
at  the  incident  bar/specimen  interface,  and  partially  reflected  at  the  specimen/ 
transmitter  bar  interface.  Strain  Gage  B  (Figure  2)  on  the  transmitter  bar 
records  the  portion  of  the  wave  that  has  transmitted  the  specimen  (oj),  while 
Strain  Gage  A  on  the  incident  bar  records  that  portion  of  the  wave  reflected 
at  the  incident  bar/specimen  interface  (oj^).  From  these  strain  gage  measure¬ 
ments,  the  stress  and  strain  in  the  specimen,  which  is  sandwiched  between  the 
two  pressure  bars,  can  be  computed  as  a  function  of  time  using  simple  wave 
mechanics. 

From  one-dimensional  theory  of  wave  analysis,  the  particle  velocity 
(V)  and  stress  (a)  in  the  bars  are  related  through  the  impedance  (i.e,,  pC^ 
where  p  is  the  mass  density  of  the  bars  and  is  the  rod  wave  velocity): 
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Figure  2.  Schematic  of  Split  Hopkinson  Pressure  Bar 
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(1) 


V  = 


The  net  particle  velocity  of  the  incident  bar-specimen  interface  is  the  result 
of  both  the  incident  and  reflected  wave. 


V 


1 


(2) 


Note  that  the  sign  of  is  opposite  that  of  Oj  and  Oj.  The  particle  velocity 
of  the  transmitter  bar-specimen  interface  is 


V 


2 


(3) 


By  averaging  the  particle  velocities  at  the  specimen-bar  interfaces,  the 
average  stress  and  strain-rate  in  the  specimen  can  be  determined  from 


(Oj  +  +  Oj)A^ 

‘^AVG  " 


(4) 


and 


^AVG 


(5) 


Where  A-j  and  1^2  3*"®  areas  of  the  pressure  bars  and  specimen  respectively, 
and  Lq  is  the  initial  length  of  the  specimen.  The  average  specimen  strain  is 
computed  by  integrating  Equation  (5). 

2.  Numerical  Analysis 

High  strain-rate  mechanical  testing  is  complicated  by  the  effects  of 

3  -1 

stress  wave  propagation.  At  strain  rates  above  10  s  it  is  difficult  to 
achieve  uniform  loading  conditions  over  the  gage  length  of  a  standard  tensile 
specimen  because  there  may  be  insufficient  time  to  dampen  the  often  complex 
stress  waves  generated  during  the  test.  The  complex  geometry  associated  with 
grips,  specimen  design,  screw  threads,  etc.,  makes  analysis  of  stress  wave 
propagation  in  such  a  test  virtually  intractable. 

The  SHPB  has  evolved  into  a  useful  high-rate  test  apparatus  because 
the  stress  waves  generated  in  long  cylinders  are  relatively  simple  and  are 
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capable  of  precise  analysis.  In  addition,  specimen  dimensions  have  been  reduced 
significantly  to  minimize  delays  associated  with  stress  wave  propagation. 
However,  several  significant  shortcomings  are  associated  with  SHPB  experiments. 
First,  it  is  not  possible  to  accurately  determine  the  stress  condition  in  the 
specimen  at  failure  from  the  available  data,  and  second,  it  is  frequently  not 
possible  to  ascertain  the  mode  of  failure  in  the  specimen.  A  comprehensive 
numerical  analysis  was  conducted  on  various  SHPB  experiments  to  gain  some  in¬ 
sight  into  those  phenomena. 

Because  of  the  dynamic  nature  of  the  impulse  loading  associated  with 
the  SHPB  experiments,  and  the  highly  nonlinear  behavior  of  the  concrete  test 
specimens,  the  finite  element  method  (FEM)  of  analysis  was  employed  in  the 
research  effort  through  the  implementation  of  the  ADINA  (Reference  17)  computer 
programs.  The  analyses  were  conducted  on  two-dimensional,  axisymmetric  models 
comprised  of  nine-node  isoparametric  finite  elements  (Reference  18).  Both 
linear  and  nonlinear  analyses  were  performed. 

C.  SCOPE 

Two  different  types  of  SHPB  experiments  were  simulated  in  the  FEM 
analyses:  (1)  Splitting-tensile  tests  and  (2)  direct  tension  tests.  In  the 
splitting  tension  analyses,  the  FEM  model  included  the  2-inch  (51  mm)  diameter 
cylinder  and  a  10-inch  (254  mm)  length  of  the  transmitter  bar.  Three  different 
load  cases  were  investigated.  In  the  direct  tension  study,  both  a  square  notch 
and  saddle  notch  specimen  were  analyzed.  The  FEM  model  included  the  entire 
lengths  of  the  incident  and  transmitter  bars,  in  addition  to  the  2-inch  (51  mm) 
diameter  specimen.  One  load  case  for  each  specimen  type  was  investigated. 
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SECTION  II 


SPLiniNG-TENSILE  TESTS 


A.  INTRODUCTION 

The  splitting- tensile  test  has  recently  established  itself  as  a  measure  of 
the  tensile  strength  of  concrete  (Reference  19).  In  a  standard  (static) 
splitting-tensile  test,  a  concrete  cylinder  of  diameter  D  and  length  L  is 
placed  with  its  longitudinal  axis  horizontal  between  the  platens  of  a  testing 
machine  as  illustrated  in  Figure  3.  The  load  is  increased  until  failure  by 
splitting  along  the  vertical  diameter  takes  place.  For  any  compressive  load  P 
on  the  cylinder,  an  element  near  the  center  on  the  vertical  diameter  of  the 
cylinder  is  subjected  to  a  vertical  compressive  stress  of 


=  IL 

z  ~  ttLD  ^z(D-z) 


1] 


(6) 


and  a  horizontal  tensile  stress  of 


,  ,  2P 
^y  ■  ttLD 


(7) 


The  horizontal  stress,  a^,  on  a  section  through  the  vertical  diameter  is 
shown  in  Figure  4  (Reference  20).  The  stress  is  expressed  in  terms  of  2P/itLD. 
It  is  observed  that  a  high  horizontal  compressive  stress  exists  in  the  vicinity 
of  the  loads.  However,  since  this  is  accompanied  by  a  vertical  compressive 
stress  of  comparable  magnitude,  a  state  of  biaxial  stress  is  produced.  There¬ 
fore,  failure  in  compression  does  not  occur. 

To  investigate  the  effects  of  strain  rate  on  the  tensile  strength  of 
concrete,  splitting-tensile  tests  of  plain  concrete  specimens  were  conducted 
on  a  Split  Hopkinson  Pressure  Bar,  which  is  illustrated  in  Figure  1.  The 
specimen  arrangement  for  the  splitting- tensile  tests  is  illustrated  in 
Figure  5.  The  cylindrical  specimens  were  2  inches  (51  mm)  in  diameter  and 
2  inches  (51  mm)  in  length.  The  static,  linear  material  properties  for  the 
specimens  were  calculated  as  follows:  The  static  compressive  strength,  f'^^  = 
7000  psi  (48.3  MPa);  the  static  tensile  strength  f^^  =  560  psi  (3.86  MPa); 
Young's  Modulus,  E  =  5.5  x  10®  psi  (37,9  GPa);  and  mass  density,  p  =  .0006747 
Ib-sec^/in^  (2.4  kg/m^). 
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Figure  3.  Illustration  of  Static  Split  Cylinder  Test 


Tension  Compression 


Horizontal  Stress  (Ty  xS!=5. ) 

Figure  4.  Distribution  of  Horizontal  Stress  (o  )  on  the 

Vertical  Diameter  for  a  Static  Spl it'^Cyl inder  Test. 
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PLAN  VIEW 


SIDE  VIEW 


Figure  5.  Splitting- Tensile  Test  Arrangement  in  the  SHPB 
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SHPB  tests  were  conducted  for  three  different  loading  conditions.  The 
stress  vs.  time  histories  for  these  cases  are  Illustrated  In  Figures  6,  7,  and 
8.  It  Is  assumed  that  the  peak  dynamic  tensile  stress,  f^^,  of  the  split 
cylinder  is  proportional  to  the  peak  transmitted  stress,  Oj,  through  the  closed 
form  solution  of  Neville  (Reference  20): 


-  IL 

td  "  ttLD 


In  which 


(8) 


(9) 


where  L  Is  the  specimen  length,  D  is  the  specimen  diameter,  and  R  is  the  radius 
of  the  SHPB. 

•  • 

Additionally,  the  loading  rate,  a,  and  the  strain  rate,  e,  in  the  specimen 
can  be  estimated  from  the  expressions 


and 


e  = 


a 

E 


(11) 


where  x  is  the  time  lag  between  the  start  of  the  transmitted  stress  wave  and 
the  maximum  transmitted  stress  (which  is  determined  from  the  stress  histories 
presented  in  Figures  6,  7,  and  8).  Table  1  summarizes  the  results  obtained 
from  the  SHPB  tests. 

To  ascertain  the  stress  condition  in  the  material  specimens  at  failure, 
a  comprehensive  finite  element  method  (FEM)  study  was  conducted  on  the  SHPB 
experiments.  Both  linear  and  nonlinear  analyses  were  performed.  From  the 
results  of  the  numerical  analyses,  the  dynamic  states  of  stress  occurring  in 
the  splitting-tensile  specimens  prior  to  failure  were  researched,  as  well  as 
the  modes  of  failure. 
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-Tensile  Test  Data  Trace  for  Load  Case 


Figure  7.  Spl i tti ng- Tcnsi Ic  Test  Data  Trace  for  Load  Case 
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Figure  8.  Splitting- Tensi le  Test  Data  Trace  for  Load  Case 


TABLE  1.  SUMMARY  OF  SHPB  RESULTS. 


Load 

Case 

No. 

Incident 

Stress 

Transmitted 

Stress 

Dynamic 

Tensile 

Stress 

Loading 

Rate 

Strain 

Rate 

Experimental 
Dynami c 
Increase 
Factor 

Oj  (psi) 

o-j-  (psi) 

ftd  (psi) 

(^  (^) 
'sec' 

(sec'b 

!td 

^ts 

(MPa) 

(MPa) 

(MPa) 

(— ) 
'sec' 

1 

-8740 

-1270 

635 

9.6  X  10® 

1.7 

1.14 

-(60.22) 

-(8.76) 

(4.38) 

(66.25) 

2 

-3040 

1520 

2.1  X  10^ 

3.8 

2.71 

-(72.76) 

-(20.97) 

(10.48) 

(144.83) 

3 

-4080 

2040 

4.2  X  10^ 

7.7 

3.64 

-(264.27) 

-(28.14) 

(14.07) 

(289.66) 

B.  LINEAR  ANALYSES 

1 .  The  FEM  Model 

An  illustration  of  the  FEM  model  employed  in  the  study  is  depicted  in 
Figure  9.  The  cylinder  is  comprised  of  1200  eight-node,  two-dimensional  finite 
elements.  To  avoid  the  development  of  artificial  reflected  stresses  at  the 
interface  between  the  cylinder  and  the  transmitter  bar  by  the  imposition  of  a 
rigid  boundary,  a  10-inch  (254  mm)  length  of  the  transmitter  bar  was  incor¬ 
porated  into  the  FEM  model.  This  portion  of  the  model  is  comprised  of  200 
eight-node,  two-dimensional  finite  elements.  The  loading  on  the  split 
cylinder  is  applied  at  the  boundary  of  the  cylinder  and  the  incident  bar, 
which  is  designated  as  the  top  of  the  cylinder. 

2.  Static  Analysis 

To  calibrate  the  FEM  model  and  verify  its  accuracy,  static  analyses 
of  the  cylinder  portion  of  the  model  were  conducted  for  the  three  load  cases 
summarized  in  Table  1.  The  statically  applied  load  was  taken  as  the  product 
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9.  FEM  Model  of  Cylinder  and  Portion  of 
Transmitter  Ban 


of  the  incident  stress  and  the  cross-sectional  area  of  the  incident  bar.  The 
numerical  results  were  compared  with  a  closed-form  analytical  solution 
(Reference  20). 

The  distribution  of  the  horizontal  and  vertical  stresses  along  the 
vertical  diameter,  obtained  from  the  FEM  analysis  for  Load  Case  1,  are  pre¬ 
sented  in  Figures  10  and  11,  respectively.  Similar  plots  for  the  horizontal 
and  vertical  stress  distributions  along  the  horizontal  diameter  are  presented 
in  Figures  12  and  13.  The  numerical  results  for  the  horizontal  and  vertical 
stresses  occurring  at  the  center  of  the  cylinder  for  all  three  load  cases  are 
presented  in  Table  2.  Excellent  correlation  with  the  closed- form  solutions  is 
noted. 


TABLE  2.  STATIC  STRESSES  AT  CENTER  OF  CYLINDER. 


Load 

Case 

No. 

Vertical  Stress, 

a,  (psi) 

^  (MPa) 

Horizontal  Stress, 

Oy  (psi) 
^  (MPa) 

FEii  Analysis 

Eq.  (5) 

FEM  Analysis 

Eq.  (6) 

1 

-1861 

-1905 

613 

635 

-(12.83) 

-(12.14) 

(4.23) 

(4.38) 

2 

-4456 

-4560 

1468 

1520 

-(30.73) 

-(31.45) 

(10.12) 

10.48) 

3 

-5980 

-6120 

1970 

2040 

-(4.124) 

-(42.21) 

(13.39) 

(14.07) 

3.  Dynamic  Analysis 

Dynamic  analyses  were  conducted  on  the  SHPB  splitting- tensile  speci¬ 
mens  described  in  the  previous  section.  The  loading  conditions  for  the  dynamic 
analyses  were  determined  from  the  SHPB  data  curves  of  the  incident,  reflected, 
and  transmitted  strain  gage  traces  presented  in  Figures  6,  7,  and  8.  These 
curves  were  corrected  for  dispersion  and  phase  change.  The  stress  on  the 
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Figure  11.  Vertical  Stress  Distribution  along  the 
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incident  face,  STRESS  1,  the  stress  on  the  transmitted  face,  STRESS  2,  and  the 
average  of  these  two  stresses,  AVE  STRESS,  are  shown  in  Figures  14,  15,  and  16, 
for  Load  Cases  1,  2,  and  3,  respectively. 

Experience  with  SHPB  experiments  (References  10,  11)  suggests  that 
the  STRESS  2  curve  is  indicative  of  the  load  transmitted  to  the  specimen.  In 
the  numerical  analyses,  the  load  functions  were  simulated  with  the  ramp  load¬ 
ing  depicted  in  Figure  17.  The  rise  time,  t^,  the  stress  level,  P^,  and  the 
time  of  duration,  t^,  for  the  three  load  cases  are  summarized  in  Table  3. 

TABLE  3.  PARAMETERS  FOR  RAMP  LOAD  FUNCTION 


Load  Case 

No. 

Rise  Time 

t^  (ysec) 

Stress  Level 

Pq  (psi) 

(MPa) 

Time  of  Duration 

t^  (psec) 

1 

66 

-1270 

100 

-(8.76) 

2 

72 

-3040 

-(20.97) 

100 

3 

48 

-4080 

100 

-(28.14) 

For  problems  in  which  an  elastic  body  is  subjected  to  a  short- 
duration  impulse  loading,  the  propagation  of  stress/strain  waves  through  the 
body  must  be  considered  in  formulating  the  solution.  Modal  analyses  generally 
do  not  yield  cost-effective,  accurate  results  for  wave  propagation  problems, 
therefore,  a  direct  numerical  integration  procedure  must  be  utilized. 

In  the  present  study,  the  Newmark  method  of  implicit  time  integration 
with  a  consistent  mass  formulation  is  employed.  The  dynamic  equilib>^ium  equa¬ 
tions  for  the  system  are  expressed  as 

[M]{U(t)}  +  [C]{LI(t)}  +  [K]{U(t)}  =  {R(t)}  (12) 
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Figure  14.  Corrected  SHPB  Da 


igure  15.  Corrected  SHPB  Data  for  Load  Case 


24 


Figure  16.  Corrected  SHPB  Data  for  Load  Case 


Time 


Figure  17.  Typical  Ramp  Loading  Function  used  in 
FEM  Analysis. 
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when  [M],  [C],  [K]  are  the  mass,  damping,  and  stiffness  matrices;  {R}  is  the 

•  •• 

external  load  vector;  and  {U},  {U},  and  {U}  are  the  displacement,  velocity, 
and  acceleration  vectors  of  the  finite  element  assemblage.  In  the  Newmark 
method,  the  following  assumptions  are  used  (Reference  21): 

^{U}  +  ^  (^{U}  +  (13) 

t+At{^}  t^-j  ^  ^  (t^"j  ^ 

In  addition  to  Equations  (13)  and  (14),  for  solution  of  the  displace¬ 
ments,  velocities,  and  accelerations  at  time  t  +  At,  the  equilibrium  equations 
(Equation  (12))  at  time  t  +  At  are  also  considered: 

[M]  +  [C]  +  [K]  (15) 


Solving  from  Equation  (14)  for  in  terms  of  ,  and  then  substitut¬ 
ing  for  into  Equation  (13),  the  equations  for  and  are 

obtained,  each  in  terms  of  the  unknown  displacements  only.  These  two 

relations  for  and  are  substituted  into  Equation  (15)  to  solve 

for  after  which,  using  Equation  (13)  and  Equation  (14),  and 

can  also  be  calculated. 

The  time  step  selected  for  temporal  integration  in  a  wave  propagation 
problem  is  critical  to  the  accuracy  and  stability  cf  the  solution.  Since  the 
Newmark  method  is  unconditionally  stable,  selection  of  the  time  step  can  be 
based  entirely  upon  accuracy.  In  a  wave  propagation  problem,  the  maximum  time 
step  is  related  to  wave  speed  in  the  material  and  element  size.  The  maximum 
time  step  is  selected  so  that  the  stress  wave  propagates  the  distance  between 
element  integration  points  within  that  time  increment.  The  maximum  time  step 
is  defined  by 


(At) 


max 


(16) 


where  is  the  length  of  an  element  in  the  direction  of  wave  propagation,  and 
c  is  the  velocity  of  wave  propagation,  given  by 
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E 

P 


(17) 


c  = 


It  has  been  determined  from  experience  (Reference  22)  that  a  time 

step  of 


yields  accurate  results.  In  the  present  study  a  time  step  of  At  =  50  nano¬ 
seconds  was  used  for  all  dynamic  analyses. 

Time  histories  for  the  horizontal  stress,  ay,^^'for4.oad  Case  1  at  five 
locations  along  the  vertical  diameter  are  illustrated  in  Figures  18,  19,  20,  21, 
and  22  for  z  equal  to  0.125D,  0.30D,  0.50D,  0.715D,  and  0.915D,  respectively. 
Similar  time  histories  for  Load  Case  2  are  presented  in  Figures  23  through  27, 
and  for  Load  Case  3  are  presented  in  Figures  28  through  32.  These  time  his¬ 
tories  indicate  that  the  maximum  horizontal  stress,  (^^y)^^^’ 
vicinity  of  the  center  of  the  cylinder  (@  z  =  0.50D).  The  values  for  (cry)^^,^^* 
for  the  three  load  cases,  are  690  psi  (4.76  MPa),  1650  psi  (11.38  MPa),  and 
2450  spi  (16.90  MPa),  as  illustrated  in  Figures  20,  25,  and  30,  respectively. 
This  observation  suggests  that,  under  dynamic  loading,  the  cylinder  would  ini¬ 
tiate  cracking  somewhere  in  its  interior  between  the  loading  face  and  middepth. 
The  crack  would  then  propagate  along  the  vertical  diameter  toward  the  outer 
boundaries  of  the  cylinder  and  eventually  perpetuate  failure. 

Time  histories  for  the  vertical  stress,  o^,  in  the  transmitter  bar 
are  illustrated  in  Figures  33,  34,  and  35  for  Load  Cases  1,  2,  and  3,  respec¬ 
tively.  The  computed  maximum  vertical  stress,  (cf2^max’  three  load 

cases  are  1540  psi  (10.62  MPa),  3400  psi  (23.45  MPa),  and  5750  psi  (34.66  MPa), 
respectively.  The  experimentally  measured  transmitted  stresses  are  illustrated 
in  Figures  6,  7,  and  8  for  Load  Cases  1,2,  and  3,  respectively.  The  corres¬ 
ponding  values  for  (a,)„^  for  the  three  load  cases  are  1300  psi  (8.97  MPa), 

3200  psi  (22.07  MPa),  and  4000  psi  (27.59  MPa),  respectively.  For  each  case 
the  measured  transmitted  stress  is  higher  than  the  calculated  stress.  However, 
this  is  to  be  expected  in  a  linear  analysis  where  no  cracking  of  the  specimens 
is  considered.  Results  for  the  same  computed  stress  based  upon  a  nonlinear 
analysis  are  in  closer  agreement  with  the  experimentally  measured  transmitted 
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Figure  18.  Time  History  for  Horizontal  Stress  at  z  =  .125D,  Load  Case 
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Figure  20.  Time  History  for  H orizontal  Stress  at  z  =  0.50D,  Load  Case 


T  I  I  I  I  r 

q  Q  o  q  Q  o 

u>  't  «  CM  ^  d 

ledWI^^O 


31 


T — 

I 

r 

in 

O 

in 

o 

C4 

P 

N 

in 

r“ 

o’ 

d 

32 


hszo 


Figure  24.  Time  History  for  Vertical  Stress  at  z  =  0.30D,  Load  Case 
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Figure  29.  Time  History  for  Vertical  Stress  at  z  =  0.30D,  Load  Case 
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Figure  32.  Time  History  for  Vertical  Stress  at  z  =  0.915D,  Load  Case 


Vertical  Stress  in  the  Transmitter  Bar,  Load  Case 


Figure  34.  Time  History  for  the  Vertical  Stress  in  the  T ransmitter  B ar.  Load  Case 
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stresses.  Moreover,  the  transmitted  stresses  depend  on  the  assumed  loading 
stress  applied  to  the  top  of  the  cylinder. 

Profiles  for  the  horizontal  stress,  o^,  at  selected  times  (see 
Table  4)  along  the  vertical  diameter  for  Load  Case  1  are  illustrated  in 
Figures  36  and  37.  Profiles  for  the  horizontal  stress  along  the  horizontal 
diameter  are  illustrated  in  Figures  38  and  39.  Similar  profiles  for  the  hori¬ 
zontal  stress  along  the  vertical  diameter  for  Load  Case  2  are  illustrated  in 
Figures  40  and  41,  and  in  Figures  42  and  43  for  Load  Case  3.  Profiles  for  the 
horizontal  stress  along  the  horizontal  diameter  are  illustrated  in  Figures  44 
and  45  for  Load  Case  2  and  in  Figures  46  and  47  for  Load  Case  3.  Examination 
of  these  profiles  reveals  the  close  resemblance  between  the  dynamic  stress 
profiles  and  the  corresponding  static  stress  profiles  presented  in  Figures  10 
and  12. 

The  representative  times  for  the  stress  profiles  are  presented  in 

Table  4.  The  numerical  results  for  the  maximum  dynamic  vertical,  and 

z  nia  X 

horizontal,  {o  )  ,  stresses  occurring  at  the  center  of  the  cylinder  (z  =  0.50D) 

y  max 

are  summarized  in  Table  5.  Also  presented  in  Table  5  are  the  dynamic  increase 
factors  (DIF)  for  each  load  case 


TABLE  4.  SELECTED  TIMES  FOR  STRESS  PROFILES  (LINEAR  ANALYSIS) 


Designation 

Time  (usee) 

Designation 

Time  (usee) 

^1 

3.14 

^6 

18.85 

^2 

6.28 

^7 

25.14 

^3 

9.43 

^8 

50.17 

^4 

12.57 

tg 

78.86 

^5 

15.71 
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Horizontal  Stress  along  the  Vertical  Diameter 


Figure  42.  Profiles  f 
Load  Case 


Profiles  for  the  Horizontal  Stress  along  the  Vertical  Diameter 
Load  Case  3*  • 
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the  Horizontal  Stress  along  the  Horizontal  Diameter 
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Figure  47.  Profiles  for  the  Horizontal  Stress  along  the  Horizontal  Diameter 


TABLE  5.  DYNAMIC  STRESSES  AT  CENTER  OF  CYLINDER 


Dynamic  Impact  Factors  (DIF) 

Load 

Case 

No. 

Vertical 

Horizontal 

^^2  ^dynamic  ^“^y^  dynamic 

^^y^ dynamic 

a2(psi) 

(MPa) 

cry(psi) 

(MPa) 

^'^z^static  ^^y^static 

1 

-2050 

690 

1.10 

1.09 

1.23 

-(14.14) 

(4.76) 

2 

-  4750 

1650 

1.07 

1.09 

2.95 

-(32.76) 

(11.38) 

3 

-7500 

2450 

1.25 

1.20 

4.38 

-(51.72) 

(16.90) 

C.  NONLINEAR  ANALYSIS 
1 .  Background 

A  nonlinear  material  analysis  was  conducted  to  ascertain  the  failure 
pattern  for  the  dynamic  splitting-tensile  tests.  The  concrete  material  model 
employed  in  the  nonlinear  analysis  was  a  hypoelastic  model  based  on  a  uniaxial 
stress-strain  relation  (Figure  48)  that  was  generalized  to  take  biaxial  and 
triaxial  conditions  into  account.  The  model  employed  three  basic  features  to 
describe  the  material  behavior:  (1)  a  nonlinear  stress-strain  relation  in¬ 
cluding  strain  softening  to  allow  for  weakening  of  the  material  under  increasing 
compressive  stresses;  (2)  a  failure  envelope  that  defines  cracking  in  tension 
and  crushing  in  compression;  and  (3)  a  strategy  to  model  postcracking  and 
crushing  behavior  of  the  material. 

An  appropriate  failure  envelope  must  be  employed  to  establish  the 
uniaxial  stress-strain  law  accounting  for  multi-axial  stress  conditions. 

Since  failure  of  the  split  cylinder  is  tension-dominated,  the  tension  failure 
envelope  depicted  in  Figure  49  was  used  in  the  concrete  model.  To  identify 
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Figure  48.  Uniaxial  Stress-Strain  Relation  used 
in  Concrete  Mode 
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Figure  49.  Three  nirnensicnal  Tensile  Failure 
Envelope  of  Concrete  Model 
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whether  the  material  has  failed,  the  principal  stresses  are  used  to  locate 
the  current  stress  state  in  the  failure  envelope.  The  tensile  strength  of  the 
material  in  a  principal  direction  does  not  change  with  the  introduction  of 
tensile  stresses  in  other  principal  directions.  However,  the  compressive 
stresses  in  the  other  principal  directions  alter  the  tensile  strength.  The 
pertinent  material  parameters  for  the  failure  envelope  and  the  uniaxial  stress 
strain  relation  are  summarized  in  Fable  6. 

TABLE  6,  CCNCRtTE  MODEL  PARAMETERS. 


Parameter 


Specified  Value 


Eq,  initial  tang. .it  modulus,  psi  (GPa) 

c^,  uniaxial  cut-off  tensile  strength,  psi  (MPa) 

a  ,  uniaxial  maximum  compressive  stress,  psi  (MPa) 

uniaxial  ultimate  compressive  si.»ess,  psi  (MPa) 

a',  uniaxial  compressive  failure  stress  under 
multiaxial  conditions,  psi  (MPa) 

e  ,  compressive  strain  at  ' 

e^,  uniaxial  ultimate  compressive  strain 

0^1  >  O  r,,  a  principal  stresses  in  airections 
pi  Pc  pj  2,  respectively 

al,  uniaxial  cut-off  tensile  stress  under 
multiaxial  conditions 


5,500,000.0  (37.93) 

550.0  (3.86) 

7,000.0  (48.28) 
7,000.0  (48.28) 

9,100.0  (62.76) 
0.0022 
0.005 

(Figure  49) 

(Figure  49) 


Because  of  the  complexity  of  the  material  description  used  for  the 
FEM  model,  an  appropriate  strategy  for  solving  the  nonlinear  finite  element 
equations  was  selected,  specifically,  the  Newton-Raphson  iteration  scheme.  In 
the  Newton-Raphson  formulation,  the  equilibrium  conditions  at  time  t  +  At  are 
satisfied  by  succes.'^ive  approximations  of  the  form  (Reference  23) 

[K]^"'{AU}''  =  {R}  -  {n^'\  (19) 

in  which  [K]^~^  is  the  tangent  stiffness  matrix  at  the  iterat'  i  -  1  and 
time  t  +  At;  {AU}^  is  the  i^^  correction  to  the  current  displacement  vector; 


62 


{R}  is  the  externally  applied  load  vector;  is  the  force  vector  that 

corresponds  to  the  current  element  stresses.  The  displacement  increment  cor¬ 
rection  is  used  to  obtain  the  next  displacement  approximation  (Reference  21) 

{U}^  =  +  {AU}\  (20) 

Equations  (19)  and  (20)  constitute  the  Newton-Raphson  solution  of  the  equili- 

O 

brium  equations  subjected  to  the  initial  conditions  [K(t  +  AT)]  =  [K(t)], 
{F(t  +  At)}  =  {^(t)},  and  {U(t  +  At)}  =  {U(t)}.  The  iteration  continues 
until  appropriate  convergence  criteria  are  satisfied. 

Using  the  Newton-Raphson  iteration,  the  governing  equilibrium  equa¬ 
tions  (neglecting  the  effects  of  a  damping  matrix  for  the  sake  of  clarity) 
presented  in  Equation  (12)  become 

[M]^+At{U}i  +  *[K]{AU}^  =  ^"^^^{R}  -  t+At^pji-l  (21) 

t+At{u}i  3  t+At(yji-l  ^  (22) 

Using  the  relations  in  Equations  (21)  and  (22),  along  with  the  assumptions 
employed  in  Equations  (13)  and  (14)  for  the  Newmark  method,  results  in 

t+At{”}i  ^  _4^(t+At(^jji-l  ^  t(|jj  ^  t("j  (23) 

At 

and  substituting  into  Equation  (21)  yields 

^[K]{AU}^  =  ^■^^^{R}  -  -  [M[(-^^'^^^{U}^'^  -  ^{U}) 

At 

-  --  ^{U}  -  ^{ii}  (24) 

where 

*[K]  =  ‘[K]  +  -T  [M].  (25) 

Ar 
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2.  Results  of  Nonlinear  Analysis  (Load  Case  1) 

Time  histories  for  the  horizontal  stress,  and  the  horizontal 
strain,  e^,  at  three  locations  along  the  vertical  diameter  are  illustrated  in 
Figures  50,  51,  and  52  for  z  equal  to  O.ID,  0.3D,  and  0.5D,  respectively. 
Similar  time  histories  for  the  vertical  stress,  a^,  and  vertical  strain,  e^, 
at  the  same  locations  are  presented  in  Figures  53,  54,  and  55.  Time  histories 
for  the  horizontal  stress  and  strain  at  three  locations  along  the  horizontal 
diameter  are  presented  in  Figures  56,  57,  and  58  for  y  equal  to  0.25D,  0.38D, 
and  0.5D,  respectively.  Similar  time  histories  for  the  vertical  scress  and 
strain  at  the  same  locations  are  presented  in  Figures  59,  60,  and  61,.  Time 
histories  for  the  vertical  stress  and  strain  in  the  transmitter  bar  are  pre¬ 
sented  in  Figure  62.  Tiie  maximum  stresses  and  strain  rates  predicted  at  each 
of  these  locations  are  summarized  in  Table  7. 


TABLE  7.  MAXIMUM  STRESSES  AND  STRAIN  RATES 
AT  SELECTED  LOCATIONS  'LOAD  CASE  1) 


Locati on 

^“^y^max 

psi 

(MPa) 

^^z^max 

psi 

(MPa) 

(sec)’^ 

• 

"z 

(sec)‘^ 

z  =  O.ID 

458.5 

(3.16) 

-8418 

(-58.1) 

4.466 

-14.9 

z  =  0.3D 

519 

(3.58) 

-3641 

(-25.12) 

3.555 

-  7.29 

z  =  0.5D 

520.4 

(3.59) 

-1995 

(-13.77) 

5.67 

-  6.625 

y  =  0.25D 

229 

(1.58) 

-1142 

(-7.68) 

NA 

-  3.J5 

y  =  0.38D 

NA 

-  568 
(-3.92) 

NA 

-  1.80 

y  =  0.50D 

NA 

NA 

NA 

NA 

Bar 

NA 

-1309 

(-9.03) 

NA 

-  1.0 

I 
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Figure  55.  Time  Histories  for  (a)  Vertical  Stress  and  (b)  Vertical  Strain 
at  z  =  0.5D,  Load  Case  1,  Nonlinear  analysis 
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Profiles  for  the  horizontal  stress,  and  the  horizontal  strain, 

Ey,  along  the  vertical  diameter  at  five  selected  time  increments  (56.1  ysec, 

56.3  lisec,  58.0  ysec,  66.0  ysec,  and  85.0  ysec)  are  presented  in  Figures  63 
through  67,  respectively.  The  initiation  of  the  first  crack  is  indicated  in 
Figure  63,  and  failure  of  the  cylinder  is  indicated  in  Figure  67.  Profiles 
for  the  horizontal  stress  and  strain  along  the  horizontal  diameter,  for  the 
same  selected  times,  are  presented  in  Figures  68  through  72.  The  propagation 
of  the  crack  to  the  center  of  the  cylinder  is  illustrated  in  Figure  70  at  time 
t  =  58  ysec.  Profiles  for  the  vertical  stress  and  strain  along  the  horizontal 
diameter  at  the  same  selected  times  are  presented  in  Figures  73  through  77. 

The  cracking  sequence  simulated  in  the  numerical  analysis,  from  ini¬ 
tiation  of  the  first  crack  until  failure,  is  illustrated  in  Figure  78.  The 
first  crack  occurs  along  the  vertical  diameter  at  a  location  approximately  0.5 
inches  (12.7  mm)  from  the  top  of  the  cylinder,  at  a  time  t  =  56.1  psec  (Figure 
78a).  At  a  time  t  =  58.0  ysec,  the  crack  is  observed  to  propagate  in  both 
directions  along  the  vertical  diameter  (Figure  78b),  past  midcylinder  in  one 
direction  and  0.2  inches  (5.1  mm)  from  the  top  in  the  other  direction.  At 
time  t  =  66  ysec,  the  crack  has  propagated  0.15  inches  (3.8  mm)  from  the  bottom 
of  the  cylinder  and  0.075  inches  (1.9  mm)  from  the  top  of  the  cylinder  (Figure 
78c).  Finally,  at  time  t  =  85  ysec,  failure  occurs  (Figure  78d).  The  crack 
has  nearly  propagated  through  the  entire  depth  of  the  cylinder,  and  flexural 
tensile  cracks  have  developed  at  either  end  of  the  cylinder  along  the  hori¬ 
zontal  diameter.  No  bifurcation  of  the  primary  vertical  crack  was  predicted  in 
this  simulation. 

3.  Res  of  Nonlinear  Analysis  (Load  Case  2) 

Timt  Histories  for  the  horizontal  stress,  a^,  and  the  horizontal 
strain,  e^,  at  three  locations  along  the  vertical  diameter  are  illustrated  in 
Figures  79,  80,  and  81  for  z  equal  to  O.ID,  0.3D,  and  0.5D,  respectively. 

Similar  time  histories  for  the  vertical  stress,  and  the  vertical  strain, 
e^,  at  the  same  locations  are  presented  in  Figures  82,  83,  and  84.  Time 
histories  for  the  horizontal  stress  and  strain  at  three  locations  along  the 
horizontal  diameter  are  presented  in  Figures  85,  86,  and  87  for  y  equal  to 
0.25D,  0.38D,  and  0.5D,  respectively.  Similar  time  histories  for  the  vertical 
stress  and  strain  at  the  same  locations  are  presented  in  Figures  88,  89,  and  90. 
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along  the  vertical  diameter.  Load  Case  1,  nonlinear  analysis, 
t  =  56.1  ysec. 
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Fiqure  64.  Profiles  for  (a)  horizontal  stress  and  (b)  horizontal  strain 
along  the  vertical  diameter.  Load  Case  1,  nonlinear  analysis, 
t  =  56.3  visec. 
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Figure  65.  Profiles  for  (a)  horizontal  stress  and  (b)  horizontal  strain 
along  the  vertical  diameter.  Load  Case  1,  nonlinear  analysis, 
t  =  58.0  jisec. 
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Figure  66.  Profiles  for  (a)  horizontal  stress  and  (b) 
along  the  vertical  diameter.  Load  Case  1, 
t  =  66,0  ysec. 
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Figure  67.  Profiles  for  (a)  horizontal  stress  and  (b)  horizontal  strain 
along  the  vertical  diameter.  Load  Case  1,  nonlinear  analysis, 
t  =  85.0  ysec. 
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Figure  68.  Profiles  for  (a)  horizontal  stress  and  (b)  horizontal  strain 
along  the  horizontal  diameter.  Load  Case  1,  nonlinear  analysi 


CBdW^AAO 


85 


Figure  69.  Profiles  for  (a)  horizontal  stress  and  (b)  horizontal  strain 

along  the  horizontal  diameter.  Load  Case  1,  nonlinear  analysis, 
t  =  56.3  ysec. 
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Figure  70.  Profiles  for  (a)  horizontal  stress  and  (b)  horizontal  strain 

along  the  horizontal  diameter.  Load  Case  1,  nonlinear  analysis, 
t  =  58.0  ysec. 


Profiles  for  (a)  horizontal  stress  and  (b)  horizontal  strain 
along  the  horizontal  diameter.  Load  Case  1,  nonlinear  analysis, 
t  =  66.0  gsec. 


ure  72.  Profiles  for  (a)  horizontal  stress  and  (b)  horizontal  strain 

along  the  horizontal  diameter.  Load  Case  1,  nonlinear  analysis, 
t  =  85.0  psec. 
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gure  73.  Profiles  for  (a)  vertical  stress  and  (b)  vertical  strain  along 

the  horizontal  diameter.  Load  Case  1,  nonlinear  analysis,  t  =  56.1  psec. 
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Figure  74.  Profiles  for  (a)  vertical  stress  and  (b)  vertical  strain  along 

the  horizontal  diameter.  Load  Case  1,  nonlinear  analysis,  t  =  56,3  ysec. 
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Figure  75.  Profiles  for  (a)  vertical  stress  and  (b)  vertical  strain  along 

the  horizontal  diameter.  Load  Case  1,  nonlinear  analysis,  t  =  58.0  ysec. 
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Figure  76.  Profiles  for  (a)  vertical  stress  and  (b)  vertical  strain  along 

the  horizontal  diameter.  Load  Case  1,  nonlinear  analysis,  t  =  66.0  jjsec. 
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Figure  77.  Profiles  for  (a)  vertical  stress  and  (b)  vertical  strain  along 

the  horizontal  diameter.  Load  Case  1,  nonlinear  analysis,  t  =  85.0  usee. 
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Figure  78.  Failure  pattern  for  splitting-tensile 

specimen.  Load  Case  1,  nonlinear  analysis. 
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Figure  79.  Time  histories  for  (a)  horizontal  stress  and  (b)  horizontal  strain 
0  z  =  0.10,  Load  Case  2,  nonlinear  analysis. 


histories  for  (a)  horizontal  stress  and  (b)  horizontal  strain 
0.3D,  Load  Case  2,  nonlinear  analysis. 


histories  for  (a)  horizontal  stress  and  (b)  horizontal  strain 
0.5D,  Load  Case  2,  nonlinear  analysis. 


histones  for  (a)  vertical  stress  and  (b)  vertical  strain 
O.ID,  Load  Case  2,  nonlinear  analysis. 


histories  for  (a)  vertical  stress  and  (b)  vertical  strain 
0.5D,  Load  Case  2,  nonlinear  analysis. 


istories  for  (a)  horizontal  stress  and  (b)  horizontal  strain 
0.25D,  Load  Case  2,  nonlinear  analysis. 


histories  for  (a)  horizontal  stress  and  (b)  horizontal  strain 
=  0.38D,  Load  Case  2,  nonlinear  analysis. 


histories  for  (a)  horizontal  stress  and  (b)  horizontal  strain 
=  0.5D,  Load  Case  2,  nonlinear  analysis. 


histories  for  (a)  vertical  stress  and  (b)  vertical  strain 
=  0.25D,  Load  Case  2,  nonlinear  analysis. 


histories  for  (a)  vertical  stress  and  (b)  vertical  strain 
0.38D,  Load  Case  2,  nonlinear  analysis. 


histones  for  (a)  vertical  stress  and  (b)  vertical  strain 
0.5D,  Load  Case  2,  nonlinear  analysis. 


Time  histories  for  the  vertical  stress  and  strain  in  the  transmitter  bar  are 
presented  In  Figure  91.  The  maximum  stresses  and  strain  rates  predicted  at 
each  of  these  locations  are  summarized  in  Table  8. 

TABLE  8.  MAXIMUM  STRESSES  AND  STRAIN  RATES 
AT  SELECTED  LOCATIONS  (LOAD  CASE  2) 


Location 

^°y^max 

psi 

(MPa) 

psi 

(MPa) 

(sec)"^ 

• 

(sec)"^ 

z  =  O.ID 

424.9 

(2.93) 

NA 

(  ) 

10.4 

-31.6 

z  =  0.3D 

461.0 

(3.18) 

-3948 

(-27.24) 

9.75 

-15.83 

z  =  0,5D 

479.8 

(3.31) 

-1840 

(-12.69) 

7.367 

-10.4 

y  =  .25D 

NA 

NA 

NA 

-  7.5 

y  =  .38D 

NA 

-  635.7 
(-4.39) 

NA 

-  4.1 

y  =  .50D 

NA 

NA 

NA 

NA 

Bar 

NA 

-2141 

(-14.77) 

NA 

-  2.08 

Profiles  for  the  horizontal  stress,  and  the  horizontal  strain, 

Ey,  along  the  vertical  diameter  at  four  selected  time  increments  (29.2  psec, 
31.2  psec,  35  psec,  and  45  psec)  are  presented  in  Figures  92  through  95,  re¬ 
spectively.  The  initiation  of  the  first  crack  is  indicated  in  Figure  92,  and 
failure  of  the  cylinder  is  indicated  in  Figure  95.  Profiles  for  the  horizontal 
stress  and  strain  along  the  horizontal  daimeter,  for  the  same  selected  times, 
are  presented  in  Figures  96  through  99.  The  propagation  of  the  crack  to  the 
center  of  the  cylinder  is  illustrated  in  Figure  98  at  time  t  =  35  psec. 
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histories  for  (a)  vertical  stress  and  (b)  vertical  strain 
transmitter  bar,  Load  Case  2.  nonlinear  analysis. 
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igure  92.  Profiles  for  (a)  horizontal  stress  and  (b)  horizontal  strain  along 

the  vertical  diameter.  Load  Case  2,  nonlinear  analysis,  t  =  29.2  ysec. 
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Figure  93.  Profiles  for  (a)  horizontal  stress  and  (b)  horizontal  strain  along 

the  vertical  diameter.  Load  Case  2,  nonlinear  analysis,  t  =  31.2  ysec. 
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the  vertical  diameter.  Load  Case  2,  nonlinear  analysis,  t  =  35  ysec. 
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lies  for  (a)  horizontal  stress  and  (b)  horizon 
vertical  diameter.  Load  Case  2,  nonlinear  anal 
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Figure  96,  Profiles  for  (a)  horizontal  stress  and  (b)  horizontal  strain  along 

the  horizontal  diameter.  Load  Case  2,  nonlinear  artalysis,  t  =  29.2  ysec. 
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Figure  97.  Profiles  for  (a)  horizontal  stress  and  (b)  horizontal  strain  along 

the  horizontal  diameter.  Load  Case  2,  nonlinear  analysis,  t  =  31.2  psec. 
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Profiles  for  (a)  horizontal  stress  and  (b)  horizontal  strain  along 
the  horizontal  diameter.  Load  Case  2,  nonlinear  analysis,  t  =  35  ysec. 
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Profiles  for  the  vertical  stress  and  strain  along  the  horizontal  diameter  at 
the  same  selected  times  are  presented  in  Figures  100  through  103. 

The  cracking  sequence  siirajlated  in  the  numerical  analysis,  from  ini¬ 
tiation  of  the  first  crack  until  failure,  is  illustrated  in  Figure  104.  The 
first  crack  occurs  along  the  vertical  diameter  at  a  location  approximately 
0.35  inches  from  the  top  of  the  cylinder,  at  a  time  t  =  29.2  ysec.  (Figure 
104a).  At  time  t  =  31.0  usee,  the  crack  is  observed  to  propagate  in  both 
directions  along  the  vertical  diameter  (Figure  104b).  At  time  t  =  35  usee, 
the  crack  has  propagated  along  the  vertical  diameter  through  the  center  of  the 
cylinder  in  one  direction,  and  has  nearly  reached  the  top  of  the  cylinder  in 
the  other  direction  (Figure  104c).  Finally,  at  time  t  =  45  usee,  failure 
occurs  (Figure  104d). 

During  the  failure  sequence  in  the  numerical  simulation,  four  dif¬ 
ferent  bifurcations  of  the  crack  pattern  are  observed  (see  Figure  104d).  The 
first  bifurcation  occurs  at  the  center  of  the  cylinder.  The  second  and  third 
bifurcations  occur  at  approximately  the  same  time;  one  at  the  crack  front  pro¬ 
pagating  toward  the  top  of  the  cylinder  at  a  distance  0.05  inches  (1.3  mm)  from 
the  top,  the  other  at  the  crack  front  propagating  toward  the  bottom  of  the 
cylinder,  at  a  distance  1.60  inches  (41  mm)  from  the  top.  The  cracks  formed  by 
the  top  bifurcation  eventually  propagate  to  the  top  surface  of  the  cylinder. 
However,  the  two  cracks  formed  by  the  bifurcation  occuring  below  the  center  of 
the  cylinder  begin  to  move  back  toward  each  other  before  bifurcating  once  again 
at  a  distance  of  0.4  inches  (10.2  mm)  from  the  bottom  of  the  cylinder.  The 
cracks  of  this  final  bifurcation  eventually  propagate  to  the  bottom  surface  of 
the  cylinder  at  which  time  failure  occurs. 

Time  histories  for  the  hoviaontal  stress,  and  horizontal  strain, 
Ey,  occurring  at  the  branch  of  the  crack  bifurcation  located  a  a  distance  of 
0.05  inches  (1.3  mm)  from  the  top  of  the  cylinder  are  presented  in  Figure  105. 
Similar  time  histories  for  the  vertical  stress,  and  vertical  strain, 
occurring  at  the  same  location  are  illustrated  in  Figure  106.  Time  histories 
for  the  horizontal  stress  and  strain  occurring  at  the  branch  of  the  crack  bifur¬ 
cation  located  0.4  inches  (10.2  mm)  from  the  bottom  of  the  cylinder  is  pre¬ 
sented  in  Figure  107.  Similar  time  histories  for  the  vertical  stress  and  strain 
at  the  same  location  are  presented  in  Figure  108. 
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Figure  100.  Profiles  for  (a)  vertical  stress  and  (b)  vertical  strain  along  the 
horizontal  diameter.  Load  Case  2, .nonlinear  analysis,  t  =  29.2  psec. 
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Figure  102.  Profiles  for  (a)  vertical  stress  and  (b)  vertical  strain  along  the 
horizontal  diameter.  Load  Case  2,  nonlinear  analysis,  t  =  35  ysec. 
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igure  103.  Profiles  for  (a)  vertical  stress  and  (b)  vertical  strain  along  the 
horizontal  diameter,  Load  Case  2,  nonlinear  analysis,  t  =  45  ysec. 


Time  histories  for  (a)  horizontal  stress  and  (b)  horizontal  strain 
at  root  of  bifurcation  at  top  of  cylinder.  Load  Case  2,  nonlinear 
analysis. 
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Time  histones  for  (a)  vertical  stress  and  (b)  vertical  strain 
at  root  of  bifurcation  at  top  of  cylinder.  Load  Case  2,  nonlinear 
analysis. 


Time  histories  for  (a)  vertical  stress  and  (b)  vertical  strain 
at  root  of  bifurcation  at  bottom  of  cylinder.  Load  Case  2,  nonlinear 
analysis. 


4.  Results  of  Nonlinear  Analysis  (Load  Case  3} 

Time  histories  for  the  horizontal  stress,  Oy,  at  three  locations 
along  the  vertical  diameter  are  Illustrated  In  Figures  109,  110,  and  111  for 
z  equal  to  O.ID,  0.3D,  and  0.50,  respectively.  Time  histories  for  the  verti¬ 
cal  stress,  at  two  locations  along  the  vertical  diameter  are  Illustrated 
in  Figures  112  and  113  for  z  equal  to  0.2D  and  0.5D,  respectively. 

Profiles  for  the  horizontal  stress,  Oy,  along  the  vertical  diameter 
at  five  selected  time  Increments  (27.1  ysec,  28.1  ysec,  28.9  ysec,  30.5  ysec, 
and  35.0  ysec)  are  presented  In  Figures  114  through  118,  respectively.  The 
Initiation  of  the  first  crack  Is  depicted  In  Figure  114,  and  failure  of  the 
cylinder  Is  Illustrated  In  Figure  118.  Profiles  for  the  horizontal  stress 
along  the  horizontal  diameter,  for  the  same  selected  times,  are  presented  In 
Figures  119  through  123.  The  propagation  of  the  crack  to  the  center  of  the 
cylinder  Is  Illustrated  In  Figure  121  at  time  t  =  29.6  ysec. 

The  cracking  sequence,  from  the  Initiation  of  the  first  crack  until 
failure.  Is  Illustrated  In  Figure  124.  The  first  crack  occurs  at  a  location 
approximately  0.4  Inches  (10.2  mm)  from  the  top  of  the  cylinder,  at  a  time 
t  =  27.0  ysec  (Figure  124a),  At  time  t  *  27.07  ysec,  the  crack  Is  observed  to 
propagate  In  either  direction  along  the  vertical  diameter  (Figure  124b).  At 
time  t  =  28.65  ysec,  the  crack  has  propagated  along  the  vertical  diameter 
through  the  center  of  the  cylinder  In  one  direction,  and  has  nearly  reached  the 
top  of  the  cylinder  In  the  other  direction  (Figure  124c).  Finally,  at  time 
t  =  35  ysec,  failure  occurs  (Figure  124d). 

During  the  failure  sequence  of  the  numerical  analysis,  three  different 
bifurcations  of  the  crack  pattern  are  observed  (see  Figure  124d).  The  first 
bifurcation  occurs  just  below  the  center  of  the  cylinder.  At  approximately  the 
same  time,  a  second  bifurcation  occurs  at  the  crack  front,  propagating  toward 
the  top  of  the  cylinder  at  an  approximate  distance  of  0.2  Inches  (5.1  mm)  from 
the  top.  The  cracks  formed  in  the  top  bifurcation  eventually  propagate  to  the 
top  surface  of  the  cylinder.  However,  the  two  cracks  formed  by  the  bifurcation 
occurring  just  below  the  center  of  the  cylinder  begin  to  move  back  toward  each 
other  before  bifurcating  once  again  at  a  distance  0.3  inches  (7.6  mm)  from  the 
bottom  of  the  cylinder.  The  cracks  of  this  final  bifurcation  eventually  propa¬ 
gate  to  the  bottom  surface  of  the  cylinder  at  which  time  failure  occurs. 
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listory  for  horizontal  stress  at  z  =  O.ID,  Load  Case 
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Figure  111.  Time  history  for  horizontal  stress  at  z  =  0.5D,  Load  Case 
nonlinear  analysis. 
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Profile  for  horizontal  stress  along  the  vertical  diameter 


Profile  for  horizontal  stress  along  the  vertical  diameter 
Load  Case  3,  nonlinear  analysis,  t  =  28.1  ysec. 
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SECTION  III 


DIRECT  TENSION  TESTS 

A.  INTRODUCTION 

The  direct  tension  test  has  seldom  been  used  to  evaluate  the  tensile 
strength  of  concrete.  This  is  because  of  the  difficulties  of  holding  the 
specimens  to  achieve  axial  tension  and  the  uncertainties  of  secondary  stresses 
induced  by  the  holding  devices.  Recently,  however,  direct  tension  tests  using 
a  Split  Hopkinson  Pressure  Bar  (SHPB)  have  been  successfully  conducted  (Refer¬ 
ence  24).  Two  types  of  tensile  specimens  were  tested,  a  square  notch  specimen 
(Figure  125a)  and  a  saddle  notch  specimen  (Figure  125b).  All  specimens  were 
cemented  to  the  ends  of  the  SHPB  with  a  nonepoxy  concrete  cement.  The  bar 
surfaces  and  specimen  surfaces  were  cleaned  in  a  manner  similar  to  that  used 
for  surface  cleaning  before  the  placement  of  foil -resistant  strain  gages. 

The  principles  of  operation  of  the  compressive  SHPB  (Figure  126)  are  de¬ 
tailed  in  Reference  25,  and  these  same  principles  apply  to  the  direct  tension 
SHPb.  The  configuration  of  the  SHPB  arrangement  employed  in  the  direct  tension 
study  reported  in  Reference  2^  is  also  illustrated  in  Figure  126.  The  tensile 
loading  mechanism  consists  oT  a  hollow  cylindrical  striker  bar  sliding  on  the 
compressive  transmitter  bar  (Bar  2)  of  the  SHPB.  The  striker  bar  impacts  a 
tup  threaded  into  the  end  of  what  becomes  the  tensile  incident  bar  (Bar  2).  A 
tensile  stress  wave  then  propagates  toward  the  specimen,  cemented  between  the 
two  bars. 

The  strain  gage  signals  for  a  square  notch  and  a  saddle  notch  specimen  ere 
presented  in  Figures  127  and  128,  respectively.  Knowing  that  the  transmitted 
and  reflected  signals  are  coincident  in  time,  then  the  peak  of  the  transmitted 
signal  is  observed  to  occur  during  the  rise  time  of  the  reflected  pulse.  More¬ 
over,  test  data  accumulated  from  a  recent  direct  tension  study  (Reference  26) 
indicates  that  failure  may  occur  in  the  rise  time  of  the  leading  pulse. 

To  ascertain  the  stress  condition  in  the  material  specimens  at  failure, 
a  comprehensive  finite  element  method  (FEM)  study  was  conducted  on  several  SHPB 
direct  tension  tests.  Both  linear  and  nonlinear  analyses  were  performed.  From 
the  results  of  the  numerical  analyses,  the  dynamic  states  of  stress  occurring 
in  the  direct  tension  specimens  before  failure,  the  modes  of  failure,  and  the 
times  of  failure  were  revealed. 
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Figure  125.  Direct  tension  specimens:  (a)  square 
notch,  (b)  saddle  notch. 
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Figure  126.  Schematic  of  Direct  Tension  SHPB. 
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Saddle  Notch  Test  Data  Trace 


B.  LINEAR  ANALYSIS 

1 .  The  FEM  Model 

To  provide  an  accurate  nurrerical  simulation  of  the  SHPB  direct  tension 
tests,  a  detailed  FEM  model  of  the  specimens  and  portions  of  the  incident  and 
transmitter  bars  was  constructed.  An  illustration  of  the  FEM  model  for  the  in¬ 
cident  bar  is  presented  in  Figure  129.  A  52-inch  (1321  mm)  segment  of  the  in¬ 
cident  bar  was  modeled  with  1594  eight-node  axisymmetric  elements.  A  similar 
representation  of  the  transmitter  bar  is  presented  in  Figure  130.  The  FEM  model 
of  the  square  notch  specimen  is  presented  in  Figure  131.  It  is  comprised  of 
348  eight-node,  axisummetric  elements.  The  FEM  model  of  the  saddle  notch  speci¬ 
men  is*  illustrated  in  Figure  132.  It  is  comprised  of  408  eight-node,  axisym¬ 
metric  elements.  For  each  analysis,  the  incident  and  transmitter  bars  were 
joined  with  the  appropriate  specimen  to  provide  a  continuous  FEM  model.  The 
longitudinal  axis  of  the  model  is  the  z-axis,  and  the  transverse  axis  is  the 
y-axis. 

2.  Calibration  of  FEM  Model 

To  verify  the  accuracy  of  the  FEM  model  in  simulating  wave  propagation, 
the  transmitter  bar  portion  of  the  model  was  subjected  to  a  simple  square-wave 
impulse.  The  intensity  of  the  input  wave  was  1000  psi  (6.9  MPa)  and  the  duration 
was  50  ysec.  Time  histories  for  the  longitudinal  stress,  0^,  at  four  locations 
along  the  length  of  the  bar  (z  =  1  inch  (25.4  mm),  z  =  12  inches  (305  mm), 
z  =  26  inches  (660  mm),  and  z  =  51  inches  (1295  mm);  where  the  origin  is 
assumed  at  the  input  end  of  the  incident  bar)  are  presented  in  Figures  133 
through  136.  The  results  predicted  by  the  numerical  ’.  -..lation  correspond 
closely  with  experimentally  recorded  stress  wave  traces. 

3.  Square  Notch  Test 

The  loading  condition  for  the  square  notch  test  was  determined  from 
the  stress  signal  presented  in  Figure  127.  The  load  function  employed  was  a 
modified  ramp  loading,  depicted  in  Figure  137.  The  peak  pressure  is  Pq  = 

5000  psi  (34.5  MPa),  the  rise  time  t^  =  47  ysec,  the  time  of  duration  for  uni¬ 
form  load  tj  =  153  sec,  and  the  total  load  duration  t^  =  200  ysec. 

Time  histories  for  the  longitudinal  stress,  a^,  at  four  locations 
along  the  longitudinal  centerline  of  the  specimen  are  presented  in  Figure  138. 
Time  histories  for  the  longitudinal  stress  at  the  notch  roots  are  presented  in 
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Figure  129.  FEM  Model  of  Incident  Bar 
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Figure  130.  FEM  Model  of  Transmitter  Bar 
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Figure  133.  Time  History  for  Longitudinal  Stress  at  z 


1016  mm  From  Spc. 
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Figure  134.  Time  History  for  Longitudinal  Stress  at 
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Figure  138.  Time  histor*  for  longitudinal  stress  along  longitudinal 
centerline,  square  notch  test. 


Figure  139.  Examination  of  these  time  histories  suggests  a  stress  concentrat¬ 
ing  factor  of  approximately  1.4  at  the  notch  roots. 

Profiles  of  the  longitudinal  stress,  along  a  transverse  section 
passing  through  the  notch  roots  are  illustrated  in  Figures  140  and  141  for 
three  selected  time  intervals  (t  =  270  ysec,  t  =  300  ysec,  and  t  =  420  ysec). 
These  profiles  provide  confirming  evidence  of  the  stress  concentration  at  the 
notch  roots. 

4.  Saddle  Notch  Test 

The  saddle  notch  specimen  was  subjected  to  the  same  loading  conditions 
as  the  square  notch  specimen.  Time  histories  for  the  longitudinal  stress, 
at  four  locations  along  the  longitudinal  centerline  of  the  specimen  are  pre¬ 
sented  in  Figure  142.  Time  histories  for  the  longitudinal  stress  at  four  loca¬ 
tions  along  the  outer  surface  of  the  specimen  are  illustrated  in  Figure  143. 
Examination  of  these  time  histories  suggests  a  stress  concentration  factor  of 
approximately  1.6  at  the  apex  of  the  notch. 

Profiles  of  the  longitudinal  stress,  a  ,  along  a  transverse  section 

A 

passing  through  the  center  of  the  specimen  are  illustrated  in  Figure  144  for 
three  selected  time  intervals  (t  =  260  ysec,  t  =  290  ysec,  and  t  =  350  ysec). 
These  profiles  provide  confirmation  of  stress  concentration  at  the  apex  of  the 
notch. 

C.  NONLINEAR  ANALYSIS 

1.  Introduction 

The  concrete  material  model  employed  in  the  nonlinear  analysis  is  a 
hypoelastic  model  based  upon  the  uniaxial  stress-strain  relation  depicted  in 
Figure  48.  The  tension  failure  envelope  illustrated  in  Figure  49  was  incor- 
proated  in  the  concrete  model.  The  pertinent  material  parameters  for  the 
failure  envelope  and  the  uniaxial  stress-strain  relation  are  summarized  in 
Table  6. 

2.  Square  Notch  Test 

Time, histories  for  the  longitudinal  stress,  a^,  and  the  longitudinal 
strain,  e^,  for  three  longitudinal  locations  along  the  exterior  surfaces  of  the 
specimen  are  illustrated  in  Figures  145,  146,  and  147,  for  z  =  0,  z  =  L/2,  and 
z  =  L,  respectively  (where  L  is  the  length  of  the  specimen  and  z  is  measured 
from  the  face  of  the  incident  bar).  Time  histories  for  and  at  the  same 
locations  along  the  longitudinal  axis  of  symmetry  are  presented  in  Figures  148, 
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Figure  139.  Tiine  histories  for  longitudinal  stress  at  the  notch  roots, 
square  notch  test. 
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Figure  140.  Profiles  of  longitudinal  stress  along  a 
transverse  plane  through  the  notch  root. 
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Figure  141.  Profiles  of  longitudinal  stress 
along  a  transverse  plane  through 
the  notch  root. 


163 


I  I  Q  u. 

O  6  Q  O  O  0U50U)^ 

<D  lO  CO  CM  ’*■  t 


t|8CqSS8J)S  l^fXV 


164 


Time  (ijsecj 

Time  histories  for  longitudinal  stress  along  the  longitudinal  centerl 
saddle  notch  test. 


“"I  I  '«l|  I  'XyV  ''W  •  1. 09 


in  Tf  CO  ca 


(BdW) 


165 


histories  for  longitudinal  stress  in  the  vicinity  of  the  notch 
e  notch  test. 
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Figure  144.  Profiles  of  the  long.itudinal  stress  along 
a  transverse  plane  through  the  notch. 
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Figure  145.  Time  histories  for  (a)  longitudinal  stress  and  (b)  longitudinal 
strain  along  the  exterior  surface,  nonlinear  analysis,  square 
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Figure  146.  Time  histories  for  (a)  longitudinal  stress  and  (b)  longitudinal 
strain  along  the  exterior  surface,  nonlinear  analysis,  square 
notch,  z  =  L/2. 
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Figure  147.  Time  histories  for  (a)  longitudinal  stress  and  (b)  longitudinal 
strain  along  the  exterior  surface,  nonlinear  analysis,  square 
notch,  z  =  L. 
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histories  for  (a)  longitudinal  stress  and  (b)  longitudinal  strain 
I  the  axis  of  symmetry,  nonlinear  analysis,  square  notch,  z  =  0. 


149,  and  150,  respectively.  The  fnaximum  stresses  and  strain  rates  predicted  at 
each  of  these  locations  are  summarized  In  Table  9. 


TABLE  9.  MAXIMUM  STRESSES  AND  STRAIN  RATES 
AT  SELECTED  LOCATIONS 
(SQUARE  NOTCH  TEST) 


Location 

(a  ) 

'  z'max 

psi 

• 

e 

z 

(sec'b 

(MPa) 

y  =  z  =  0 

546 

(3.77) 

5.35 

v=^  z=t 

y  2*  ^  r 

348 

(2.40) 

4.14 

y  =  z  =  L 

535 

(3.69) 

6.12 

y  a  0,  Z  =  0 

542 

(3.74) 

4.49 

y  =  0,  z  =  ^ 

555 

(3.83) 

3.2 

II 

N 

o 

II 

535 

(3.69) 

4.98 

Transmitter  Bar 

589 

(4.06) 

2.57 

Profiles  of  the  longitudinal  stress,  a^,  for  three  selected  times  at 
three  locations  transverse  to  the  longitudinal  axis  (z  =  0,  z  =  L/2,  and  z  =  L) 
are  Illustrated  in  Figures  151,  152,  and  153,  respectively.  Cracking  in  the 
specimen  at  the  root  of  the  notch  and  at  the  specimen  incident  bar  Interface  Is 
evidenced  In  Figures  152  and  151,  respectively. 

The  cracking  sequence  simulated  by  the  numerical  analysis,  from  ini¬ 
tiation  of  the  first  crack  until  failure,  is  illustrated  in  Figures  154,  155, 
and  156.  The  first  cracking  in  the  specimen  occurs  at  the  roots  of  the  notch 
at  time  t  =  270  ysec  (Figure  154).  This  is  consistent  with  the  stress 
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Figure  149.  lime  histories  for  (a)  longitudinal  stress  and  (b)  longitudinal  strain 
along  the  axis  of  symmetry,  nonlinear  analysis,  square  notch,  z  =  L/2. 
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histories  for  (a)  longitudinal  stress  and  (b)  longitudinal  strain 
I  the  axis  of  symmetry,  nonlinear  analysis,  square  notch,  z  =  L. 
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Figure  151.  Profiles  for  longitudinal  stress  transverse 


to  the  longitudinal  axis,  nonlinear  analysis, 
square  notch,  2*0. 
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Figure  152.  Profiles  for  longitudinal  stress  transverse 
to  the  longitudinal  axis,  nonlinear  analysis 
square  notch,  z  =  L/2. 
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Figure  154.  Cracking  in  square  notch  specimen  at  time  t  =  270.0  ysec. 
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Figure  155.  Cracking  in  square  notch  specimen  at  time  t  =  274.0  psec. 
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Figure  156.  Cracking  in  square  notch  specimen  at  time  t  =  275.0  psec. 


concentration  predictions  from  the  linear  analysis.  However,  at  time  t  = 

274  lisec,  cracks  develop  at  the  outside  surface  of  the  specimen  next  to  the 
incident  bar  at  t  =  274  ysec  (Figure  155).  Failure  occurs  along  a  transverse 
plane  at  the  end  of  the  specimen  adjacent  to  the  incident  bar  at  time  t  = 

275  psec  (Figure  156).  This  failure  pattern  is  not  evidenced  in  the  linear 
analysis. 

3.  Saddle  Notch  Test 

Time  histories  for  the  longitudinal  stress,  a^,  and  the  longitudinal 
strain,  e^,  for  three  longitudinal  locations  along  the  exterior  surface  of  the 
specimen  are  illustrated  in  Figures  157,  158,  and  159,  for  z  =  0,  z  =  L/2,  and 
z  =  L,  respectively  (where  L  is  the  length  of  the  specimen  and  z  is  measured 
from  the  face  of  the  incident  bar).  Tin«  histories  for  and  at  the  same 
locations  along  the  longitudinal  axis  of  symmetry  are  presented  in  Figures  160, 
161,  and  162,  respectively.  The  maximum  stresses  and  strain  rates  produced  at 
each  of  these  locations  are  summarized  in  Table  10. 

Profiles  of  the  longitudinal  stress,  o^,  for  three  selected  times  at 
three  locations  transverse  to  the  longitudinal  axis  (z  =  0,  z  =  L/2,  and  z  =  L) 
are  illustrated  in  Figures  163,  164,  and  165,  respectively.  Cracking  in  the 
specimen  at  the  apex  of  the  notch  and  at  the  specimen-incident  bar  interface 
is  evidenced  in  Figures  164  and  163,  respectively. 

The  cracking  sequence  simulated  by  the  numerical  analysis,  from  ini¬ 
tiation  of  the  first  crack  until  failure,  is  illustrated  in  Figures  166,  167, 
and  168.  The  first  cracking  in  the  specimen  occurs  at  a  transverse  section 
in  the  specimen,  adjacent  to  the  face  of  the  incident  bar,  at  a  time  t  =  265 
ysec  (Figure  160).  This  observation  is  contradictory  to  the  anticipated  fail¬ 
ure  location  predicted  by  the  linear  analysis.  At  time  t  =  257  ysec,  signifi¬ 
cant  crack  growth  along  this  same  plane  is  indicated  (Figure  164).  Finally, 
failure  occurs  along  this  plane  at  time  t  =  258  ysec  (Figure  168).  It  is  not 
until  time  t  =  258  ysec  that  cracks  develop  in  the  notch,  contemporaneous  with 
failure  at  another  location.  This  failure  pattern  is  not  supported  by  the 
results  of  the  linear  analysis. 
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Time  histories  for  (a)  longitudinal  stress  and  (b)  longitudinal 
strain  along  the  exterior  surface,  nonlinear  analysis,  saddle  notch 


ir  (a)  longitudinal  stress  and  (b)  longitudinal 
exterior  surface,  nonlinear  analysis,  saddle  notch 
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Figure  159.  Vime  histories  for  (a)  longitudinal  stress  and  (b)  longitudinal 
strain  along  the  exterior  surface,  nonlinear  analysis,  saddle 

notch,  z  =  L. 
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Figure  160.  Time  histories  for  (a)  longitudinal  stress  and  (b)  longitudinal 
strain  along  the  axis  of  symmetry,  nonlinear  analysis,  saddle 


JO  A' 


'|U|«J)8f1)  AAa 


|BdM|  AAo 


Figure  161.  Time  histories  for  (a)  longitudinal  stress  and  (b)  longitudinal 
strain  along  the  axis  of  symmetry,  nonlinear  analysis,  saddle 
notch,  z  =  L/2. 
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Figure  162.  Time  histories  for  (a)  longitudinal  stress  and  (b)  longitudinal 
strain  along  the  axis  of  sytmietry,  nonlinear  analysis,  saddle 
notch,  z  =  L. 
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Figure  163.  Profiles  for  longitudinal  stress  transverse  to 
the  longitudinal  axis,  nonlinear  analysis, 
saddle  notch,  z  =  0. 
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Figure  164.  Profiles  for  longitudinal  stress  transverse  to 
the  longitudinal  axis,  nonlinear  analysis, 
saddle  notch,  z  =  L/2. 
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Figure  165.  Profiles  for  longitudinal  stress  transverse  to 
the  longitudinal  axis,  nonlinear  analysis, 
saddle  notch,  z  -  L. 
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Figure  168.  Cracking  in  saddle  notch  specinien  at  time  t  =  258.0  psec. 


TABLE  10.  MAXIMUM  STRESSES  AND  STRAIN 

RATES  AT  SELECTED  LOCATIONS 
(SADDLE  NOTCH  TEST) 

Location 

psi 

• 

^Z 

(sec’b 

(MPa) 

y  =  2  =  0 

552 

(3.81) 

9.62 

y  =  D  2  =  L 

y  2*  2 

558 

(3.85) 

21.8 

y  =  z  =  L 

554 

(3.82) 

90.4 

y  =  0.  z  =  0 

557 

(3.84) 

8.88 

y  *  0,  z  =  1 

556 

(3.84) 

8.80 

y  =  0,  z  *  L 

559 

(3.86) 

9.84 

Transmitter  Bar 

538 

(3.71) 

2.36 
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SECTION  IV 


DISCUSSION  OF  RESULTS 

A.  SPLiniNG-CYLINDER  ANALYSIS 

1.  Linear  Analysis 

The  results  of  the  linear  analyses  indicate  that  the  dynamic  stress 
distribution  in  the  cylinder  behind  the  initial  stress  wave  is  identical  to 
that  exhibited  in  the  static  analysis.  This  fact  is  particularly  noticeable 
in  the  figures  illustrating  the  profiles  for  the  horizontal  stress,  along 
both  the  vertical  and  horizontal  diameters  (Figures  36  through  47).  This 
observation  suggests  that  a  dynamic  failure  would  closely  resemble  the  static 
failure.  However,  the  results  of  the  nonlinear  analyses  dispell  this  assump¬ 
tion. 

Another  interesting  observation  pertains  to  the  development  of  the 
maximum  horizontal  stress  at  the  center  of  the  cylinder.  Again,  this  is  con¬ 
sistent  with  the  results  of  the  static  analysis.  This  fact  can  be  shown  from 
the  aforementioned  profiles  of  the  horizontal  stress,  a^,  and  from  the  time 
histories  of  at  selected  locations  along  the  vertical  diameter  presented  in 
Figures  18  through  32.  Based  upon  these  results,  one  could  conclude  that  the 
cylinder  would  fail  from  the  center,  outward  to  the  exterior  boundaries  of  the 
vertical  diameter. 

2.  Nonlinear  Analysis 

The  mode  of  failure  predicted  by  the  nonlinear  analysis  differs  from 
that  suggested  by  the  results  of  the  linear  analysis.  In  all  three  load  cases, 
the  initiation  of  first  cracking  is  not  at  the  center  of  the  cylinder  (as  the 
results  of  the  linear  analysis  indicate),  but  at  an  approximate  distance  of 
0.2D  from  the  top  of  the  cylinder.  The  failure  sequence  for  Load  Cases  1,  2, 
and  3  are  illustrated  in  Figures  78,  104,  and  124,  respectively.  The  sequence 
of  failure  for  all  load  cases  is  essentially  the  same:  initiation  of  first 
cracking  at  location  0.2D  from  the  top,  and  subsequent  propagation  of  the 
cracks  in  both  directions  along  the  vertical  centerline  of  the  cylinder  toward 
the  top  and  bottom  surfaces. 

The  deviation  of  the  failure  pattern  predicted  by  the  nonlinear 
analysis  from  that  predicted  by  the  results  of  the  linear  analysis  can  be 
attributed  to  the  tension  failure  envelope  employed  in  the  analysis.  The 
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tension  failure  envelope  is  depicted  in  Figure  49.  To  identify  material  failure, 
the  principal  stresses  are  used  to  locate  the  current  stress  state  in  the 
failure  envelope.  The  tensile  strength  of  a  material  in  a  principal  direction 
does  not  change  with  the  introduction  of  tensile  stresses  in  other  principal 
directions,  however,  the  compressive  stresses  in  other  principal  directions 
alter  the  tensile  stress.  Since  the  compressive  stresses  in  the  vicinity  of 
the  top  of  the  cylinder  are  high  in  comparison  to  those  at  the  center  of  the 
cylinder,  initiation  of  first  cracking  will  be  at  that  location. 

Another  interesting  idiosyncrasy  of  the  mode  of  failure  is  revealed 
in  the  illustrations  of  the  failure  sequences  presented  in  Figures  78,  104,  and 
124.  The  failure  pattern  for  Load  Case  1  (Figure *78)  reveals  no  bifurcation 
in  the  primary  crack  pattern.  For  Load  Case  2  (Figure  104),  however,  several 
bifurcations  in  the  primary  crack  pattern  are  observed.  And  for  Load  Case  3 
(Figure  124),  an  exaggerated  bifurcation  pattern  is  observed.  It  can  be  con¬ 
cluded  that  the  presence  and  extent  of  the  bifurcations  in  the  failure  pattern 
are  related  to  the  load  rate.  Load  Case  1  represents  a  relatively  low  load 
rate,  while  Load  Cases  2  and  3  represent  successively  high  load  rates.  There¬ 
fore,  it  is  concluded  that  the  load  rate  affects  the  mode  of  failure. 

The  patterns  of  cracking  predicted  by  the  numerical  analysis  is  con¬ 
sistent  with  those  observed  in  the  SHPB  experiments  (Reference  24).  Preliminary 
results  of  high  speed  photography  (10,000  frames/second)  taken  at  AFESC  of  high 
load  rate  SHPB  tests  indicate  the  development  of  cracks  similar  to  the  pattern 
illustrated  in  Figures  104  and  124,  along  the  vertical  diameter  prior  to  the 
appearance  of  the  bifurcated  cracks  at  the  top  and  bottom  of  the  cylinder. 
Moreover,  the  SHPB  specimens  have  exhibited  evidence  of  crack  bifurcation 
occurring  just  below  the  center  of  the  cylinder,  similar  to  that  illustrated  in 
Figure  124d.  This  is  substantiated  by  the  observation  from  high-speed  photo¬ 
graphy  of  a  lens-shaped  piece  of  fractured  concrete  being  expelled  from  the 
flat  surface  of  the  specimen.  It  should  be  noted  that  all  material  fractures 
are  predicated  upon  the  failure  envelope  presented  in  Figure  49.  No  fracture 
mechanics  parameters  are  used  to  describe  the  fracture  process. 

The  results  of  the  splitting-cylinder  analyses  may  be  used  to  corre¬ 
late  the  dynamic  tensile  strength  of  the  concrete  to  load  rate  when  viewed  in 
conjunction  with  experimental  strength  versus  strain  rate  data  associated  with 
the  appropriate  loading  rates. 
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B.  DIRECT  TENSION  TESTS 

1.  Linear  Analysis 

The  results  of  the  linear  analysis  indicate  the  development  of  high 
stress  concentrations  at  the  root  of  the  notch  (1.4)  for  the  square  notch  speci¬ 
men,  and  at  the  apex  of  the  notch  (1.6)  for  the  saddle  notch  specimen.  These 
concentration  factors  are  determined  from  the  time  histories  for  the  longitudi¬ 
nal  stress,  illustrated  in  Figures  138  and  139  for  the  square  notch  test, 
and  Figures  142  and  143  for  the  saddle  notch  test. 

The  results  of  the  linear  analyses  indicate  that  first  cracking  will 
begin  at  the  locations  of  high  stress  concentration.  It  can  further  be  assumed 
that  failure  will  ultimately  occur  on  vertical  planes  passing  through  those 
points. 

2.  Nonlinear  Analysis 

The  mode  of  failure  predicted  by  the  nonlinear  analysis  differs  from 
that  suggested  by  the  results  of  the  linear  analysis.  In  the  square  notch  speci¬ 
men,  first  cracking  occurs  at  the  roots  of  the  notch  (Figure  154).  This  is  con¬ 
sistent  with  the  stress  concentration  predictions  from  the  linear  analysis. 
However,  eventual  failure  of  the  specimen  occurs  on  a  vertical  p^ane  adjacent 
to  the  face  of  the  incident  bar  (Figures  155  and  156). 

In  the  saddle  notch  specitren,  first  cracking  occurs  at  a  transverse 
section  in  the  specimen  next  to  the  incident  bar  (Figure  166).  Almost  simul¬ 
taneously,  cracks  develop  in  the  apex  of  the  notch  (Figure  167).  Eventual 
failure  is  along  the  transverse  section  adjacent  to  the  incident  bar. 

The  failure  patterns  predicteo  by  the  nonlinear  analysis  contradict 
those  suggested  by  the  linear  analysis.  There  are  two  basic  reasons  for  this 
discrepancy:  (1)  the  load  rate  is  very  high,  therefore  critical  stresses  de¬ 
velop  at  the  loaded  end  of  the  specimens  before  significant  stresses  can  develop 
in  the  vicinity  of  the  notches;  and  (2)  the  notches  are  relatively  shallow, 
therefore,  they  do  not  represent  a  critical  section  for  such  a  high  rate  of 
loading. 

In  the  experimental  procedure,  the  saddle  notch  specimen  tended  to 
fail  at  both  the  bottom  of  the  saddle  notch  and  at  the  specimen  end  next  to  the 
incident  bar  for  the  same  load  rate  used  in  the  numerical  analysis.  However, 
for  the  square  notch  specimen,  the  failure  was  usually  at  the  notch  for  the 
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loading  used  in  the  numerical  analysis.  At  higher  load  rates  (i.e.,  higher 
striker  impact  velocities)  the  square  notch  specimen  also  failed  at  both  the 
notch  and  the  incident  end. 
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SECTION  V 


CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

1.  Splitting-Tensile  Tests 

The  results  of  the  linear  analysis  indicate  that  the  dynamic  stress 
distribution  in  the  cylinder  behind  the  initial  stress  wave  is  identical  to 
that  exhibited  in  the  static  analysis.  The  crack  patterns  and  the  modes  of 
failure  predicted  by  the  nonlinear  FEM  analysis  is  consistent  with  that  ob¬ 
served  in  the  SHPB  experiments. 

In  the  nonlinear  analysis, failure  by  separation  of  the  cylinders,  is 
predicted  at  time  t  =  85  psec  (Figure  78)  for  Load  Case  1,  t  =  45  psec  (Figure 
104)  for  Load  Case  2,  and  t  =  35  psec  (Figure  124)  for  Load  Case  3.  In  Load 
Cases  2  and  3,  the  failure  occurs  before  the  maximum  load  is  reached;  that  is, 
during  the  rise  time  of  the  load  (refer  to  Table  3).  This  failure  is  predi¬ 
cated  upon  a  uniaxial  tension  cut-off  stress  equal  to  the  static  tensile 
strength  of  the  material.  However,  experimental  strength  vs.  strain  rate  data 
associated  with  this  loading  rate  indicates  that  the  concrete  tensile  strength 
may  be  three  to  four  times  as  great  as  the  static  tensile  strength.  Therefore, 
using  the  dynamic  tensile  strength  in  the  material  model  would  delay  the  time 
of  failure  and  possibly  bring  it  into  the  range  of  constant  load;  that  is,  at 
a  time  greater  than  the  rise  time.  It  is  anticipated  that  the  higher  tensile 
strength  would  not  affect  the  overall  failure  pattern. 

It  can  also  be  concluded  that  the  nature  of  the  failure  mode  is 
directly  affected  by  the  rate  of  loading.  For  a  relatively  low  load  rate,  such 
as  Load  Case  1,  the  failure  mode  manifests  itself  as  a  single  crack  propagating 
along  the  vertical  centerline  of  the  cylinder  (Figure  78).  However,  for  in¬ 
creasingly  higher  load  rates,  such  as  Load  Case  2  and  Load  Case  3,  the  mode  of 
failure  is  characterized  by  several  bifurcations  in  the  primary  crack  pattern 
(Figures  104  and  124).  The  higher  the  load  rate,  the  more  pronounced  are  the 
bifurcations. 

2.  Direct  Tension  Tests 

The  results  of  the  linear  analysis  indicate  high  stress  concentrations 
in  the  vicinity  of  the  notches.  In  the  square  notch  test,  a  stress  concentration 
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factor  of  1.4  is  indicated  at  the  roots  of  the  notch.  In  the  saddle  notch  test, 
a  stress  concentration  factor  of  1.6  is  indicated  at  the  apex  of  the  notch. 

These  results  suggest  a  failure  in  a  transverse  plane  passing  through  the  notch. 

The  nonlinear  analysis,  however,  predicts  failure  in  a  transverse 
plane  near  the  end  of  the  specimen  next  to  the  incident  bar.  The  reason  for 
this  is  that  the  load  rate  is  so  high,  that  the  tensile  limit  of  the  material 
is  reached  at  the  end  of  the  specimen  (adjacent  to  the  incident  bar)  before  any 
significant  stresses  can  develop  on  the  transverse  plane  passing  through  the 
notch. 

The  mode  of  failure  is  dominated  by  the  load  rate.  Moreover,  it  can 
be  concluded  that  the  specimen  fails  during  the  rise  time  of  the  load,  before 
multiple  reflections  can  develop  within  the  specimen.  However,  it  is  possible 
that  the  failure  mode  could  be  altered  by  considering  material  strain  rate 
effects  in  the  analysis  and/or  deeper  notches  in  the  specimen. 

B.  RECOMMENDATIONS 

The  failures  predicted  in  both  the  splitting- tension  tests  and  the  direct 
tension  tests  are  highly  sensitive  to  the  rate  of  loading.  Therefore,  it  is 
recommended  that  additional  analyses  be  conducted  at  a  wide  range  of  load  rates 
to  quantify  the  relationship  of  load  rate  to  mode  of  failure. 

It  is  also  apparent  from  the  results  of  the  analyses  that  material  strain 
rate  effects  will  delay  the  time  of  failure,  allowing  the  specimen  to  be  sub¬ 
jected  to  a  higher  load,  thus  possibly  affecting  the  failure  mode.  Therefore, 
it  is  further  recommended  that  additional  numerical  analyses  be  conducted  to 
investigate  material  strain  rate  effects  on  the  mode  of  failure. 

Finally,  the  notches  in  the  direct  tension  specimens  analyzed  in  this  study 
were  relatively  shallow.  It  is  recommended  that  specimens  with  deeper  notches 
be  analyzed,  both  experimentally  and  numerically,  to  quantify  the  effect  of 
notch  depth  on  the  mode  of  failure. 
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